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^ Abstract 

o 

We consider the inverse problem of determining an optical mask that produces a 
desired circuit pattern in photolithography. We set the problem as a shape design prob- 
1^ lem in which the unknown is a two-dimensional domain. The relationship between the 

target shape and the unknown is modeled through difFractive optics. We develop a vari- 
ational formulation that is well-posed and propose an approximation that can be shown 
04 to have convergence properties. The approximate problem can serve as a foundation to 

numerical methods. 

Oh AMS 2000 Mathematics Subject Classification Primary 49Q10. Secondary 49J45, 

49N45. 

i^H Keywords photolithograpy, shape optimization, sets of finite perimeter, F-convergence. 

L_li 1 Introduction 

Photolithography is a key process in the production of integrated circuits. It is the process 
by which circuit patterns are transferred onto sihcon wafers. A review of this manufacturing 
•/^ technology is given in [16]. The main step in photolithography is the creation of a circuit 

image on the photoresist coating which sits on the silicon layer that is to be patterned. The 
image is formed using ultra-violet (UV) light which is diffracted by a mask, and refracted by 
a system of lenses. The mask simply consists of cut-outs, and lets light through the holes. 
<^ The parts of the photoresist that are exposed to the UV light can be removed, leaving 

openings to the layer to be patterned. The next stage is etching, which removes material in 
^ the layer that is unprotected by the photoresist. Once etching is done, the photoresist can 

^ be removed, and the etched away "channels" may be filled. The entire process is illustrated 

schematically in Figure 1. 

The problem we address in this work is the inverse problem of determining what mask 
is needed in order to remove a desired shape in the photoresist. The difficulty of producing 
a desired shape comes from the fact that the UV light is diffracted at the mask. Moreover, 
the chemicals in the photoresist reacts nonlinearly to UV exposure - only portions of the 
photoresist that have been exposed to a certain level of intensity are removed in the bleaching 
process. 

The nature of the present work is analytical. Our goal is to formulate mathematically 
well-posed problems for photolithography. The methods we use to prove well-posedness are 
constructive and may serve as a foundation for a computational method. 



o 



*Dipartimento di Matematica e Informatica, Universita degli Studi di Trieste, via Valerio, 12/1, 34127 
Trieste, Italy, rondi@units.it 

^ School of Mathematics, University of Minnesota, Minneapolis, MN 55455, USA. santosaSmath . umn . edu 



1 



UV light 




mask 



optics 



-photo-resist 
layer to be developed 



developed resist 



-etched layer 



resist removed 



Figure 1: The photolithographic process. Ultraviolet light, diffracted by a mask, forms an 
image on the photoresist. The exposed portion of the photoresist is removed, leaving open- 
ings. Etching removes parts of the layer to be patterned. After etching, the photoresist is 
removed. 



Our investigation into photolithography is inspired by the work of Cobb [3] who was 
the first to approach this problem from the point of view of optimal design which utilizes 
a physically-based model. This general approach was further developed by introducing a 
level set method in |18j . A different computational approach which models the mask as a 
pixelated binary image can be found in |14j . 

The plan of the paper is as follows. In the first and preliminary section. Section [2j 
we develop the most basic model for removal of the exposed photoresist. We describe the 
inverse problem to be solved. This is followed by a discussion of the approximate problem 
whose properties we intend to investigate in this work. Section [3] contains mathematical 
preliminaries needed for our work. We introduce the basic notation and recall various results 
which will be useful for our analysis. In particular, in Subsection 3.3 we discuss the geometry 
of masks or circuits and how to measure the distance between two of them. In Section [4] 
we discuss the properties of the operator which maps the mask into the circuit. Section [5] 
provides an analysis of the variational approach to the problem of the optimization of 
the mask and we prove a convergence result for it, Theorem 5.5, in the framework of F- 
convergence. 



2 Description of the inverse problem 

This section is separated into three subsections. First, we review some basic facts about 
Fourier transforms and prove a result about approximation of a Gaussian. We follow this 
with a discussion of the optics involved and a model for photolithography. In the final 
subsection we describe the inverse problem and its approximation. 
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2.1 Fourier transform and approximation of Gaussians 

We first set some notation and describe a few preliminary results. For every a; G M^, we 
shall set x = (xi,a;2), where xi and X2 G M. For every x G M'^ and r > 0, we shall denote 
by Bj.{x) the open ball in centered at x of radius r. Usually we shall write instead of 
5^(0). We recall that, for any set C M^, we denote by xe its characteristic function, and 
for any r > 0, Br{E) = \J^^eBt{x). 

For any / G 5'(M^), the space of tempered distributions, we denote by / its Fourier 
transform, which, if / G L^(M^), may be written as 

/(O = / /(x)e-'«-dx, e G 
We recah that f{x) = {2tt)-'^ f{-x), that is, when also / G L^(]R2), 

/(^) = 7^2/ my^-^dC, X G 
If / is a radial function, that is f{x) = 0(|x|) for any x G M^, then 

/(O = 2^?^o(<A)(iei), e G 

where 

/•+00 

T-Lo{(p){s) = / r Jo(sr)(/)(r)dr, s > 0, 
Jo 

is the Hankel transform of order 0, Jq being the Bessel function of order 0, see for instance 
a- 

We denote the Gaussian distribution by G{x) = {2tt) ^e 1^' x G M^, and let us note 
that G{S,) = e~l^l ^ G M^. Moreover, ||G||^i(]k2) = 1. Furthermore if 60 denotes the Dirac 

delta centered at 0, we have 6q = 1, therefore (27r)~^i = 6q. 

For any function / defined on and any positive constant s, we denote fs{x) = 
s-^f{x/s), x G M?. We note that ||/.||ii(M2) = ||/||ii(M2) and M) = f{si), ^ G M^. 

We conclude these preliminaries with the following integrability result for the Fourier 
transform and its applications. 

Theorem 2.1 There exists an absolute constant C such that the following estimate holds 

II/IIli(ir2) < C'll/llvy2.i(R^)- 

Proof. This result is contained in Theorem A in [9] and it is based on previous analysis 
done in [13]. □ 

We recall that a more detailed analysis on conditions for which integrability of the 
Fourier transform holds may be found in [TTj. However the previous result is simple to use 
and it is enough for our purposes, in particular for proving the following lemma. 

Lemma 2.2 For any 5 > there exist a constant sq, < sq < 1, and a radial function 
f G C^{M?) such that f = 1 on B^^ and, if we call T = {27ry^f, then T G 1^2,1(^2) 

- Clival. i(R2) < 6. 
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Proof. We sketch the proof of this result. Let us consider the following cut-off function 
(j) G C°^(M) such that (j) is nonincreasing, (j) = 1 on(— oo, 0] and (p = on [1, +oo). 

We define a function T as follows 

f{x) = - 1) + (1 - (j){\x\ - l))GZ{xmx\ -b), x£ 

for suitable constants sq, < sq < 1 and b>2. We call f = (27r)~^r. 



Then lengthy but straightforward computations, with the aid of Theorem 2.1, allow us 
to prove that for some sq small enough and for some b = SQ^bo, with bo large enough, we 
have 

\\T - Gs(,||vi/i,i(K2) < 6. 
Then, let T{x) = T(x/so), x G M^, so that T = Ti/^^, or equivalently T = T^p. Therefore 

Gen IIm/I-ICIBZ-i < S. 



By a simple rescaling argument we have that T satisfies the required properties. Further- 
more, by this construction, we may choose T such that it is radially nonincreasing, T = 1 
on Bgg and it decays to zero in a suitable smooth, exponential way. □ 



2.2 A model of image formation 

We are now in the position to describe the model we shall use. The current industry standard 
for modeling the optics is based on Kirchhoff approximation. Under this approximation, the 
light source at the mask is on where the mask is open, and off otherwise (see Figure 1). 
Propagation through the lenses can be calculated using Fourier optics. It is further assumed 
that the image plane, in this case the plane of the photoresist, is at the focal distance of the 
optical system. If there were no diffraction, a perfect image of the mask would be formed 
on the image plane. Diffraction, together with partial coherence of the light source, acts to 
distort the formed image. 

The mask, which we mention consists of cut-outs, is represented as a binary function, 
i.e., it is a characteristic function of the cut-outs. Suppose that D represents the cut-outs, 
then the mask is given by 

m{x) = xd{x). 

The image is the light intensity on the image plane. This is given by |12j 

(2.1) I{x)= / m{i)K{x - i)J{i-ri)K{x -r])m{r])didri, x G M^. 

In the above expression the kernel K{-) is called the coherent point spread function and 
describes the optical system. For an optical system with a circular aperture, once the 
wavenumber of the light used. A; > 0, has been chosen, the kernel depends on a single 
parameter called the Numerical Aperture, NA. Notice that the wavelength is A = 27r/A;. Let 
us recall that the so-called Jinc function is defined as 

Jinc(x) = ^^iM, ^eM2, 
2tt\x\ 
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where Ji is the Bessel function of order 1. We notice that in the Fourier space, see for 
instance P, page 14], 

jhrc(c) = xBi(o, ^£^2. 

If we denote by s = (fcNA)^^, then the kernel is usually modeled as follows 

K{x) = Jincs(x) = V , ' x£W, 

27r \x\ 

therefore 

If NA goes to +00, that is s — )• 0^, then K converges pointwise to 1, thus K approximates 
in a suitable sense the Dirac delta. 

For technical reasons, we shall consider a slightly different coherent point spread function 
K. Let us fix a positive constant 5, to be chosen later. We shall replace the characteristic 
function XBi, the Fourier t rans form of the Jinc function, with the function T[sq^), ^ G M?, 



with T and sq as in Lemma 2.2, Therefore T(sq-) is a radial function that it is still identically 
equal to 1 on i^i, it is still compactly supported, it is nonincreasing with respect to the radial 
variable and it decays to zero in a smooth, exponential way. Its Fourier transform is Tg^ 
and we shall assume that 

(2.2) K{x) = {Ts,)s{x)=Tss,{x), xeR\ 

where again s = {kNA)~^. Also in this model, if NA goes to +00, that is s — t- 0"*", then K 
converges pointwise to 1, thus K approximates in a suitable sense the Dirac delta. 

The function J(-) is called the mutual intensity function. If the illumination is fully 
coherent, J = 1. In practice, illumination is never fully coherent and is parametrized by a 
coherency coefficient a. A typical model for J is 

(2.3) Jix) = 2^^^^^=Tr,-]mcikaNA\x\), x G M^. 



Thus, 



(27r)2'^^^^ 7r(A;aNA)2^^^''^^^^^' 



that, as cr — )• 0"*", converges, in a suitable sense, to the Dirac delta. Therefore full coherence 
is achieved for cr —t- 0^. In fact, if o" — t- O"*" , J converges to 1 uniformly on any compact subset 



of M^. The equation (2.1) is often referred to as the Hopkins areal intensity representation. 
As it will become apparent from the analysis developed in the paper, the value of s is related 
to the scale of details that the manufacturing of the mask allows, thus in turn to the scale 
of details of the desired circuit. Therefore, we typically consider fcNA » 1, that is s ^ 1, 
and ka^A < 1. 

2.3 The inverse problem and its approximation 

The photoresist material responds to the intensity of the image. When intensity at the 
photoresist goes over a certain threshold, it is then considered exposed and can be removed. 
Therefore, the exposed pattern, given a mask m(x), is 

(2.4) O = {x G : I{x) > h}, 
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where h is the exposure threshold. Clearly, Q depends on the mask function m{x), which 
we recall is given by the characteristic function of D representing the cut-outs, that is 
Q = i}{D). In photolithography, we have a desired exposed pattern which we wish to 
achieve. The inverse problem is to find a mask that achieves this desired exposed pattern. 
Mathematically, this cannot, in general, be done. Therefore, the inverse problem must be 
posed as an optimal design problem. 

Suppose the desired pattern is given by Oq- We pose the minimization problem 



(2.5) 



min d{Q{D),no) 



The distance function d{-, •) will be discussed in detail below. The admissible set A is our 
search space, and needs to be defined carefully as well. 



Instead of solving (2.5), we pose a variational problem for a function u (instead of the 
mask D). We will show below that this problem is well-posed and that as the approximation 



parameter is set to zero, we recover the solution of (2.5) under a perimeter penalization. 



Instead of dealing with the characteristic function xd{x) which represents the mask, 
we will work with a phase-field function u which takes on values of and 1 with smooth 



transitions. Thus, the intensity in (2.4) is calculated with u instead of m = xd in (2.1 ), so / 



is a function of u. At this point, we will not be precise about the space of functions to which 
u belongs. To force u to take on values of mostly and 1, we introduce the Mordica-Mortola 
energy 



1 



W{u) + e 



where W{t) = 9t'^{t - if is a double-well potential. We will regularize the problem of min- 
imizing the distance between the target pattern and the exposed region by this energy. 

Then we relax the hard threshold in defining the exposed region Q in (2.4). Let (j){t) 
be a C°° nondecreasing approximate Heaviside function with values (j){t < —1/2) = and 
(^(t > 1/2) = 1. The function 

'I{u) - 

7] 



$^(n) 



will be 1 where the intensity / > h + rj/2. A sigmoidal threshold function is employed in 
the computational work in |14j. 



Now we consider the distance function between O and $7o in (2.5). Let 



(2.6) 



(i = d(fi,Qo) 



/ 



Xn-Xn,\ + \Pm-P{^o) 



where is the characteristic function of the set Vt and P{^) is the perimeter of the region 
O. To approximate this distance function, we replace it by 



+ 



\V{%{u))\-P{no] 



The characteristic function of is replaced by the smooth threshold function while its 
perimeter is replaced by the TV-norm of the function. 
The approximate problem we shall solve is 



min . 
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The remainder of the paper is an analytical study of this minimization problem. We will 
show that it is well-posed, and that in the limit e — )• O"*", we recover the solution of the 
original problem (2.5) under a perimeter penalization. 



3 Mathematical preliminaries 

By Ti^ we denote the 1-dimensional Hausdorff measure and by we denote the 2-dimensio- 
nal Lebesgue measure. We recall that, if 7 C M-^ is a smooth curve, then restricted to 7 
coincides with its arclength. For any Borel E cM? we denote \E\ = C'^{E). 

Let P be a bounded open set contained in M^, with boundary dV. We say that V has a 
Lipschitz boundary if for every x = (xi, X2) € dD there exist a Lipschitz function : M — t- M 
and a positive constant r such that for any y £ Br{x) we have, up to a rigid transformation, 

y = (yi,y2) e if and only if y2 < ^{yi). 

We note that T> has a finite number of connected components, whereas &D is formed by a 
finite number of rectifiable Jordan curves, therefore n^{dV) = length(9P) < +00. 

We recall some basic notation and properties of functions of bounded variation and sets 
of finite perimeter. For a more comprehensive treatment of these subjects see, for instance, 

miziii. 

Given a bounded open set V C M^, we denote by BV{T>) the Banach space oi functions of 
hounded variation. We recall that u G BV{T>) if and only if u S L}iT>) and its distributional 
derivative Du is a bounded vector measure. We endow BV{T>) with the standard norm 
as follows. Given u G BV(T>), we denote by \Du\ the total variation of its distributional 
derivative and we set ||ti||By(x)) = ll^llLip) + l-^^K^)- We shall call P{u,'D) = {DuKV). 
We recall that whenever u G W^'^{'D), then u € BV{V) and \Du\{V) = |Vti|, therefore 

IkllsyCD) = IkllLicc) + l|V'u||ii(x)) = 

We say that a sequence of BV{V) functions {ufej^i weakly* converges in BV{'D) to 
u G BV(T>) if and only if converges to u in and Duh weakly* converges to Du in 

D, that is 

(3.1) lim / vdDuh = / vdDu for any v E Co{'D). 

By Proposition 3.13 in [Ij, we have that if a sequence of BV{T)) functions {u^}^^-^ is 
bounded in BV(T>) and converges to u in L^(V), then u G BV{V) and Uh converges to u 
weakly* in BV{V). 

We say that a sequence of BV{T>) functions {uh}^^i strictly converges in BV{T>) to 
u G BV(T>) if and only if converges to u in L^{T>) and {DuhKl)) converges to |Z)ii|(P). 
Indeed, 

dst{u,v)= / \u - v\ + \\Du\{V) - \Dv\{V)\ 
Jv 

is a distance on BV{'D) inducing the strict convergence. We also note that strict convergence 
implies weak* convergence. 

Let V he a bounded open set with Lipschitz boundary. A sequence of BV{'D) functions 
Wh}hLi such that sup/j ||n/j||^y(-p) < +00 admits a subsequence converging weakly* in 
BV{T>) to a function u G BV{'D), see for instance Theorem 3.23 in [1]. As a corollary, we 
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infer that for any C > the set {u £ BV{'D) : ||ti||_By(X') < C} is a compact subset of 
L\V). 

For any fixed constant R > 0, with a shght abuse of notation, we shall identify L^{Bji) 
with the set {u G L-^(M^) : u = Q a.e. outside Br\. 

Let ii^ be a bounded Borel set contained in C M^. We shall denote by xe its charac- 
teristic function. We notice that E is compactly contained in -B/j+i) which we shall denote 
hy E (S: Brj^i. We say that is a set of finite perimeter if xe belongs to BV{B nj^i) 
and we call the number P{E) = IDxeKBr-^-i) its perimeter. Analogously, for any u G 
L^{Bji) n BV{Bji^i), we shall denote P{u, Bfj^i) = \Du\{Br^i). Obviously, if u = xe, 
then P{u,Br+i) = P{E). 

Let us further remark that the intersection of two sets of finite perimeter is still a set 
of finite perimeter. Moreover, whenever E is open and 7i^{dE) is finite, then E is a set of 
finite perimeter, see for instance [71 Section 5.11, Theorem 1]. Therefore a bounded open 
set D with Lipschitz boundary is a set of finite perimeter and its perimeter P{T>) coincides 
with n^{dV). 

3.1 F-convergence approximation of the perimeter functional 

Let us introduce the following, slightly different, version of a F-convergence result due to 
Modica and Mortola, [llj. We shall follow the notation and proofs contained in [2j. We 
begin by setting some notation. For the definition and properties of F-convergence we refer 
to J5j. 

For any bounded open set V C M?, with a slight abuse of notation, we identify Wq'^{V), 
1 < p < +00, with the subset of W^'^{M?) functions u such that u restricted to V belongs 
to Wq'P{V) and u is equal to almost everywhere outside P. Let us assume that for some 
positive constant R we have P C Br. We recall that any function in -L^(I') is extended 
to zero outside V and the same procedure is used for L}{Bji). Therefore, with this slight 
abuse of notation, L^{V) C L^{Bji). Throughout the paper, for any p, I < p < +oo, we 
shall denote its conjugate exponent by p' , that is -|- = 1. 

Theorem 3.1 Let V C Br C be a bounded open set with Lipschitz boundary. Let us 
also assume that D is convex. 

Let 1 < p < +00 and VF : M — t- [0, -|-oo) be a continuous function such that W{t) = if 
and only if t G {0, 1}. Let Cp = {J^ {W{s))^/p' ds)'^ . 

For any e > we define the functional Pe : L^{E?) — )• [0, -|-oo] as follows 

' ^ [ W{u) + ^'^ f \Vu\P ifueW''P{V), 

(3.2) Pe{u) = i Jv P Jv 

+ oo otherwise. 

Let P : L1(M2) [0, +oo] be such that 

(3.3) P{u) 

+ oo otherwis 
Then P = F-lim£^o+ Pe with respect to the L^{M.'^) norm. 



' P{u,Br+i) ifue BV{Br+i), u G {0,1} a.e., 
and u = a.e. outside T>, 
+ oo otherwise. 
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Remark 3.2 We observe that P{u) = P{E) if n = xe where E is a. set of finite perimeter 
contained in D and P{u) = +00 otherwise. 

Furthermore, we note that the result does not change if in the definition of we set 
Pe{u) = +00 whenever u does not satisfy the constraint 

(3.4) < li < 1 a.e. in V. 

Proof. We sketch the proof following that of Theorem 4.13 in [2j . In fact, the only diff'erence 
with respect to that theorem is that we assume P convex and that we take Wq'^{V) instead 
of W'^'P{V) in the definition of P^. 

By Proposition 4.3 in [2], we obtain that P{u) < r-liminf£^o+ Pe{u) for any u G L^(M?). 
In order to obtain the P-limsup inequality, we follow the procedure described in Section 4.2 
of 12]. It would be enough to construct Ai C L^{M?) such that the following two conditions 
are satisfied. First, we require that, for any u G L^(M?) such that P{u) < +00, there exists 
a sequence {uj}'jLi such that Uj E Ai, for any j G N, uj — )• u in L^{M?) as j — )• 00, and 
P{u) = limj P{uj). Second, for any u G Ai, r-limsup£^o+ -fsl^) ^ 

We choose M. = {u = xe '■ E V, E of class C°°}. The second property follows by 
Proposition 4.10 in [2]. As far as the first property is concerned, this can be obtained by 
following the proof of Theorem 1.24 in [8]. That theorem states that any bounded set of 
finite perimeter E can be approximated by a sequence of C°° sets {Ej}'jLi such that, as 
j — )• 00, J^2 \XEj — Xe\ and P{Ej) — )• P{E). If we assume that E C V, and that V is 
convex, by choosing in the proof of Theorem 1.24 in [8J a value of t satisfying 1/2 < t < 1, 
we obtain that the sets Ej are also compactly contained in D, for any j £N. □ 

Also the following result, due to Modica, pO], will be useful. 

Proposition 3.3 For any C > 0, let us take 1 < p < +00 and any e > 0, and let us define 

Ac = {ue L^(]R2) : < u < 1 a.e. and Ps{u) < C}. 
Then Ac is precompact in L^{M.'^). 

Proof. We repeat, for the reader's convenience, the arguments developed in |10j . Clearly 
Ac is a bounded subset of L^{'D). Let {n^j^i be a sequence in Ac. We need to prove 
that there exists a subsequence converging in L^{T)). For any t, < t < 1, let (f){t) = 
jl{W{s)Y/p'ds. For any n G N, we define Vn = (piun) and we observe that < w„ < 0(1) 
almost everywhere. Therefore, the functions u„, n G N, are uniformly bounded in L°°{T)) 
and, consequently, in L^(2?). Furthermore, since (/> is a function, with bounded norm, 
then Dvn = 4>'{un)Dun = W^/P {un)Dun- Therefore, 

/ \DVn\ = [ \W'/P'{Un)\\DUn\<Pe{Un)/Cp. 

Jv Jv 

We infer that there exists a subsequence {vnk}'kLi converging, as k ^ 00, to a function vq 
in L}{T>) and almost everywhere. Let if) be the inverse function of <j) and let uq = iP{vq). We 
observe that ijj is bounded and uniformly continuous on [0,0(1)], hence we conclude that, 
as A; — )• 00, Uuf. converges to uq in L^{V). □ 

Remark 3.4 With the same proof, we can show the following. Let us consider any family 
{ue}o<£<so such that, for some positive constant C and for any e, < e < eoj we have 
< Ue < 1 almost everywhere and Peiu^) < C. Then {ue}o<£<eo is precompact in L^(M^). 
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3.2 Convolutions 

We recall that, for any two functions / and g defined on M^, we define the convolution of / 
and f * g, as follows 

{I*9){x)= f{x -y)g{y)dy = f{y)g{x - y)dy, x G 

whenever this is well-defined. 

The following classical properties of convolutions will be used. First convolution is com- 
mutative. Second, as a consequence of Young inequality we have the following result about 
integrability and regularity of convolutions. 

Proposition 3.5 Let I < r, p, q < +00 be such that 1 + ^ = ^ + |> '^'^d let n = 0,1,2, .. .. 

Letl<q< +00, let / G L«(M2) and let g G ^^"'^(M^). Then h = f*ge W^^'iR^) and 
there exists a constant C, depending on n, p, q and r only, such that 

Let q = +00 and let f G L°°(R^), with compact support. If g ^ W"''^{M.'^), then h = 
f * g & 1/F"'°°(M^) and there exists a constant C , depending on n only, such that 

ll^llw".°°(M2) < C'II/IIl°°(m2)||5||w.i(]r2). 

If f e L^(M2) and g G L'^{M.'^), then h = f *g e L°°{M.'^) and it holds ||/i||Loo(R2) < 
||/||^oo(i[j2) ||(7||^i(]g2) . Furthermore, if g is uniformly continuous and cOg denotes its modulus 
of continuity, then h is also uniformly continuous and 

i^h < ||/||Li(]R2)a;(,. 

Finally, let f G L^{M.^) and let g G ^'"'{M.^), for some a,0<a<l. Then h G C"'°^(M2) 
and there exists a constant C, depending on n and a only, such that 

II^I|C".°(R2) < C'||/||Ll(R2)||5'||cn,a(]R2). 

3.3 The geometry of masks and circuits 

In this subsection we investigate the following two questions, namely what arc reasonable 
assumptions on the geometry of the mask D and how to measure the distance between 
the constructed circuit O and the desired one Qq. We begin with the following definition. 
During this subsection, in most cases proofs will be omitted and left to the reader. 

For given positive constants r and L, we say that a bounded open set C is 
Lipschitz or C*^'^ with constants r and L if for every x G dVt there exists a Lipschitz 
function </? : M ^ M, with Lipschitz constant bounded by L, such that for any y G Br{x), 
and up to a rigid transformation, 

(3.5) y = (y^^y^) en if and Only if y2 < ^{yi)- 

Without loss of generality, we may always assume that x = (0, 0) and i^(0) = 0. We shall 
always denote by ei and ^2 the vectors of the canonical bases. Clearly the orientation of the 
canonical bases may vary depending on x G dfl. 
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We shall also use the following notation. There exist positive constants 6i < 1/2, 82 < Si 
and mi < 1, all of them depending on L only, such that the following holds. For any x £ 
and for any S > 0, let Ms{x) = {y : \yi\ < 5r, 2/2 = fiui)} and Ns{x) = {y : \yi\ < 
6ir, ^p{yi) — 5r < y2 < ^(yi) + 6r}. Then we assume that, for any 5, < 6 < 82, the 
following properties hold. First, Ns{x) C B^/2i^) (hence Ms^{x) C B^i2{x) as well). Clearly 
Ns{x) is contained in Bsr{dQ), and we assume that Ns{x) contains Bmi6r{^5i/2{x)) and 
that for any y £ {y : \yi\ < 5ir/2, 1/2 = V'ivi) ± M) V ^ Bm^Srid^)- 

For any integer k = 1,2,..., any a, < a < 1, and any positive constants r and L, we 
say that a bounded open set C is C^'°' with constants r and L if for every x G dO, 
there exists a C^'" function 99 : M — )• M, with C^'"' norm bounded by L, such that for any 



y £ Br{x), and up to a rigid transformation, (3.5) holds. Without loss of generality, we may 
always assume that x = (0, 0) and '^{0) = 0. 

Let us fix three positive constants r, L and R. Let A^'^{r, L, R) be the class of all 
bounded open sets, contained in Bf( C M?, which are Lipschitz with constants r and L. For 
any integer k = 1,2, . . . and any a, < a < 1, we denote with A^'"'{r, L, R) the class of all 
bounded open sets, contained in Br C M^, which are C'^''^ with constants r and L. 

Since we shall identify open sets D with their characteristic functions xd, if ^ = 
A^'^{r, L, R), (or A = A'''°^{r,L,R), respectively) then, with a slight abuse of notation, 
A will also denote the subset of functions u € L^{Bfj) such that u = xd for some D £ A. 
Moreover, we shall denote 

A = {u£ L^{Br) : < u < 1 a.e. in Br} 

and, for any 7 > 0, 

(3.6) A^ = {u £ A : \\u — Xd\\l'^{Br) ^ 7 ^oi some D £ A}. 

Let us assume that and O2 belong to A^'^{r, L, R). There are several ways to define the 
distance between these two sets. We shall describe four of them and study their relationships. 
We let 

(3.7) di = di{ni,n2) = dH{ni,Q2); 

(3.8) di = dii^i,n2) = dH{dni,dn2); 

(3.9) d2 = d2{fii,^2) = \^iA^2\ = WxQi - Xnih^BR+iY, 

(3.10) d3 = d3ini,n2) = d2 + \P{ni)-Pin2)\=dstixn„Xn2)- 

Here dn denotes the Hausdorff distance, whereas we recall that P{i^) denotes the perime- 
ter of in Br^i and dgt is the distance inducing strict convergence in BV{Brj^i). First of 
all, we observe that all of these are distances. We now investigate their relationships. 

We begin with the first two, di and di, and we notice that 

(3.11) if di < r/4, then di < di. 

There exists a constant c, < c < 1, depending on L only, such that 

di > cminjr, di}. 

Therefore, 

if di < r, then di < Cdi, 
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where C = 1/c. Furthermore, if di < (c/2)r, then di must be less than or equal to r, so 

if di < (c/2)r, then di < Cdi. 

Moreover, we can find a constant ci, < ci < 1, depending on L only, such that 

if di < cir, then di < di < Cdi. 

We conclude that we can find a constant ci, < ci < 1, and a constant C > 1, both 
depending on L only, such that 

(3.12) if either di < cir or di < cir, then di < di < Cdi. 
Since di and di are bounded by 2R, we also have 

2^ 27? 

(3.13) if both di > cir and di > cir, then di < di and di < di. 

cir cir 



We finally observe that the estimates (3.12) and (3.13) are essentially optimal. 

Before comparing di (or di) with d2 and d^, let us make the following remark on the 
lengths of dili and dO,2- If ^ is an open set which is Lipschitz with constants r and L, then 
for any integer n > 0, we have 



(3.14) 



n\dQ n (5(„+i)A^nr)) < C{L)r{n + 1). 



Here, a simple computation shows that we may choose C(L) = ASViTL^. 
Therefore, if we assume that Q C Bji and R > lOr, we may conclude that 



(3.15) 



P{n) < Ci{L)R^/r, 



1 ^ii)V(L). 



where Ci(L) — 2 v g 

Moreover, there exist two constants C2, < C2 < ci, and Ci > 0, depending on L only, 
such that we have 



(3.16) 
Since 



\Bd{d^)\ < Ci\ength{d^)d for any d < C2r. 
if di < C2r, then ^2 < min{|Bj^(5Sl2)|, (SJ^s)!}, 



we obtain that 

(3.17) if di < C2r, then d2 < Ci min{length((90i), length(5r22)}di- 



If di > C2r, then d2 < ttR"^ < ^^^di. By (3.15), we may conclude that 



(3.18) 



d2 < 



CoR"" 



di. 



Here C2 d epends on L only. Moreover, up to changing the constants C2, Ci and C2, (3.17) 
and (3.18) still hold if we replace di with di. 
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On the other hand, there exists a constant C3, < C3 < vr, depending on L only, such 
that 

d2 > C3min{r^,df}. 

We infer that either if di < r or 

if d2 < (c3/2)r2, then di < C^d!^'^ , where C3 = l/cg^^ If 
^2 ^ (c3/2)r^, then d\ < 2R < ^^Rd2 or, better, di < ^^"^^ Rd}^"^ . Summarizing, we have 

(3.19) if d2 < (c3/2)r2, then di < Czd^^ 
and, finally, 

(3.20) if d2 > (c3/2)r2, then di < '^^'^^ R^j'^ . 

r 

Clearly, up to suitably changing the constants C3 and C3, the last two estimates still hold 
if we replace di with di. We also remark that, as before, the estimates relating di, di and 
d2 are essentially optimal. 

We have obtained that di, di and ^2 are topologically equivalent distances. About ^2 
and (i3, obviously ^2 < ^3, however the two distances are not topologically equivalent. In 
fact we can find and ilj, i G N, open sets belonging to J\P'^{r,L,R), such that d2{^,^i) 
goes to zero as z — )• 00, whereas (i3(0,ilj) > c > for any i G N. Therefore ^3 induces a 
strictly finer topology than the one induced by d2 

An assumption that the mask is a bounded open set which is Lipschitz with given 
constants r and L is reasonable from the manufacturing point of view as well as from the 
mathematical point of view, by the following compactness result. 

Proposition 3.6 The set ^°'^(r, L, R) {respectively ^^'°(r, L, /?), A; = 1, 2, . . ., < a < 1) 

is compact with respect to the distance di. 

We remark that the same result holds with respect to the distances di and ^2- Further- 
more, we obtain as a corollary that the set is closed with respect to the norm, for 
any 7 > 0. 

The previous example shows that compactness fails with respect to the distance ^3, at 
least for the Lipschitz case. On the other hand, if Q.i and ^2 belong to ^^'°(r, L, i?), with 
< a < 1, then, following Lemma 2.1 in [15], we can show that 

(3.21) |P(f)i)-P(f]2)| <C4(di(J7i,172))"/('°+'\ 

where C4 depends on r, L, R and a only. We may conclude that in the C'^'" case. A; = 1, 2, . . ., 
< a < 1, ^3 is topologically equivalent to the other three distances and that Proposition |3.6| 
holds also with respect to the distance ^3. 

It is worthwhile to observe that, under some circumstances, the estimate ( |3.21[ ) can 
be extended to the piecewise C^'" case. For example, typically we may assume that the 
desired circuit belongs to ^"^'^(r, L, R). Moreover, we assume that the boundary of VLq is 
composed by a finite number of closed segments = 1, . . . , n, which are pairwise internally 
disjoint and whose lengths are greater than or equal to 2r. Therefore, is actually a 
piecewise C^'" open set. We shall show in Section |4] that, under suitable assumptions on 
the mask D, the corresponding constructed circuit belongs to ^^'"(ri, Li, i?), for some 
suitable positive constants ri < r, Li > L, R > R and a, < a < 1. Then we can find 
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positive constants C4, < C4 < 1, C5 and Cq, depending on ri, Li, i? and a only, such that 
if di{Qo, 0) < c^ri, then we can subdivide d^l into smooth curves Ji, i = 1 . . . ,n, which are 
pairwise internaUy disjoint, such that for any i = 1, . . . , n we have 

dH{J^,I^) <CrJi{no,^) 

and 

length(/i) - 2C5diino,^) < length(Ji) < length(/,) + C6{diino,n)f/^^''+^\ 
Therefore, 

-2nC5dii^o,n) < p{n) - p{no) < nC6{di{no,n))''/^^''-^^\ 

By these reasonings it might seem that we may choose to measure the distance between 
the desired circuit f^o and the reconstructed one Q by using any of these distances. However, 
there are several reasons to prefer the distance ^3, which we actually choose. In fact, it is 
easier to compute than di and di, it can be extended in a natural way from characteristic 
functions to any BV function by using dgt, and should provide a better approximation of 
the desired circuit than ^2, which seems to be too weak for this purpose. 



3.4 Convolutions of characteristic functions and Gaussian distributions 

We recaU that G{x) = (27r)~ie-l^'l'/^ x £ R^, and let us note that G{C) = e'l^l'/^ ^ G M^. 
Moreover, ||G||ii(iR2) = 1. For any positive constant s we denote by Gs{x) = s~'^G{x/s), 

X £ note that \\Gs\\Li{R2) = 1 and Gs{C) = G{sC), C G I^^- 

Let D be a bounded open set which is Lipschitz with constants Rq and L and let xd 
be its characteristic function. We investigate how xd is perturbed if we convolute it with 
G. We call v = xd * G, that is 



v{x) = I xd{x - y)G{y)dy = I XD{y)G{x - y)dy, xGM. 
We recall that we shall use the positive constants 5i, 82 and mi, and the sets M^^ and Ns 



introduced at the beginning of Subsection 3.3 



Proposition 3.7 Under the previous notation and assumptions, let us fix 5 , Q < 6 < (52/4. 
Then there exist constants Rq > 1, h, < h < 1/24, and ai > 0, depending on L and 6 
only, such that the following estimates hold. For any x G M^, 

(3.22) ifh<v{x)<l-h, then x GBmiSRoidD), 
and for any x G dD, 

(3.23) if ye Ns{x), then Vv{y) ■ (-62) > ai. 
Proof. If x ^ BmiSRoidD), then we have 

v{x) < e-"^?^'«o/2^ iix^D, 

and 

v{x) > l-e-"^?'^'^o/2^ ifxGL>. 
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Consequently, provided h = e ^o/^ < 1/24, we may conclude that (3.22) holds. 

Let us take x G dD and y G Ns{x). Then, denoting by v the exterior unit normal vector 
to D, 

yv{y) ■ (-62) = / G{y- z)v{z) ■ e2dV}{z). 

JdD 

Therefore, 

Vv{y)-{-e2)= [ _ G{y - z)u{z) ■ e2dn\z)+ 
JdDnB2SRoiy) 

[ _ G{y - z)i^{z) ■ e2dn\z)+ 

[ Giy-z)uiz)-e2dn\z) = A + B + G. 

JdD\Bj,^/2iy) 

Since -B^g/2(y) is contained in Bji^^{x), for any z £ dD n B^^i2{y)i we have v[z) ■ 62 > 
ci > where ci is a constant depending on L only. Moreover, the length of dD n B2SR0 (y) 
is also bounded from below by C2SRQ, C2 > depending on L only. Therefore, we obtain 
that A > ciC2dRoe-^^''^o and B>0. 



For what concerns the term C, with the help of (3.14), we can find a constant Ci, 
depending on L only, such that, for any i?o > 1, we have 

\C\ < C7ii?oe-^o/8. 

Therefore, we can find Ro > 1, depending on L and 6 only, such that h = e~™i^^^o/2 < 
1/24, and 2Gie-^o/» < ciC2Se-^^' ^0 . We set ai = {l/2)ciC2SRoe-^^'^o and the proof is 
concluded. □ 

Remark 3.8 Without loss of generality, we may choose Rq such that it also satisfies 

(3.24) l|VG||,.i(M2\B«^^/,) < (l/12)ai. 

In the sequel, we shall fix 5 = 62/^ and Rq as the corresponding constant in Propo- 



sition 



3.7 such that (3.24) holds. We note that, in this case, 5 and Rq depend on L only. 
We shall also fix a constant R > 10i?o- We recall that, with a slight abuse of notation, we 
identify L^{Bji) with the set of real valued L^(]R^) functions that are equal to zero almost 



everywhere outside Br. The same proof of Proposition 3.7 allows us to prove this corollary. 



Corollary 3.9 For any s,0<s<l,letr = sRq and let D be a bounded open set which is 
Lipschitz with constants r and L. Let Vg = Xd * Gs- Then, for any x £M? , 

if h < Vs{x) < 1 — h, then x £ B^iSridD), 

and for any x G dD, 

if y £ Ns{x), then Vvs{y) ■ (-62) > aiRo/r = ai/s. 
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We conclude this part with the fohowing perturbation argument. Let us consider a 
function ip such that either iIj £ C^(M?) D W^'^{M?) or ^ G W'^'^{M?) and that, for some 
6>0, 



Let G = G + tp. Then the following result holds. 

Corollary 3.10 Let us assume that 6 < min{/i, ai/2}. 

For any s, < s < 1, let r = sRq and let D be a bounded open set which is Lipschitz 
with constants r and L. Let Vs = Xd * Gs- Then, for any x G M^, 

(3.25) if2h < Vs{x) <l-2h, then x G Bm,5r(.dD), 
and for any x G dD, 

(3.26) ifyeNsix), then Vvs{y) ■ (-62) > aiRo/2r = ai/{2s). 



Proof. It follows immediately from the previous corollary and Proposition 3.5 We first 
notice that in either cases xd * Gs ^ C^(M?). Moreover we have 

llxD *Gs-XD* Gs\\l^(m:^) < ||xdIIl-'(r2)IIV's||li < IIV'IIli(r2) 

and 

||V(XD *Gs-XD* Gs)||z,«=(R2) = WXD * {ViPs)\\l^{R2) < 

IIXD||L°o(R2)||V^/;s||ii < ||VV'||l1(r2)/s. 

Thus the conclusion follows. □ 

4 Relationship between a mask and its image intensity 

In this section we study the relationship between a function representing a mask (not nec- 
essarily a characteristic function of a domain) and its associated image intensity. We recall 



the notation used. We fix (5 = (^2/4 and Rq as the corresponding constant in Proposition 3.7 



such that (3.24) holds. We note that, in this case, 6 and Rq depend on L only. We shall also 
fix a constant R > WRq. We recall that, with a slight abuse of notation, we identify L^{Bfi) 
with the set of real valued L^(M'^) functions that are equal to zero almost everywhere outside 
Br. We recah that A = {u £ L^{Br) : < n < 1 a.e. in Br}. 
Fixed 5 > 0, we assume that ip G W'^'^{R'^) and that 



< 6. 



We denote G = G + i/j and, for any s, < s < 1, we define the operator Ts ■ L^{Br) — 
14^2,1(^2^ follows 

Ts{u) = u* Gs, for any u G L^{Br). 

The point spread function we use, T, can be described in general by the function G. There- 
fore a study of properties of convolutions with G will be useful. 
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Wc remark that the following continuity properties of the operator Tg hold. For any p, 
1 <p < +00, and any u G L}{Br), we have, for an absolute constant C, 

(4-1) ||'7^('")||lp(M2) < ||G||Ll(]R2)||u||ip(K2), 

(4.2) ||V7;(u)||ip(]R2) < (C/s)||VG||ii(]K2)||n||ip(]R2), 

(4.3) p2r.(n)|Up(K2) < (C/s2)||d2g|Ui(u2)||u|Up(k2). 

Let J G C°(M2). For any u G L\Br), we define [/ G L1(M4) as follows 
U{x,y) = u{x)u{y)J{x — y), for any x,y E M^. 
Then, for any s, < s < 1, we define Hg G VF^'^(M^) in the following way 

Hs{x,y) = Gs{x)Gs{y), for any x,y e M^. 

Therefore, for any p, 1 < p < +00, and any u G L^{Bji), we have, for an absolute constant 
C, 

(4.4) * -f^s||LP(R4) < ||G|||l(]g2)|| J||L<x.(B2^)||n|||p(jj2), 

(4.5) \\V{U *Hs)\\lp{r4) < (C/s)||G||ii(]K2)||VG||ii(]K2)||J||ioc(B2fl)||?^|lL(M2), 

(4.6) \\D'^{U * i/s)||LP(R4) < (C/s2)||G||^2.1(li2)|| J||loc(S2«)||^^||L(R2). 

Let us fix p > 4 and let a = 1 -4/p, < a < 1. Then, if « G A we have U*Hs G C^'"(M'^) 
and, for some absolute constant C depending on p, 

\\U * Hs\\c'^,<^{r4.) < C\\U * Hs\\w^•p{R■i)■ 
We define Vj^s ■ A C^'°'{M.'^) and Vi^s ■ A C^'"(M2) as follows. For any u G ^ 

Vj,s{u){x)= [ [ u{OGs{x-OJ{C-v)Gs{x-vHv)d^dr], xeR^ 

and 

ViA^) = {Ts{u)f. 
We notice that the two definitions are consistent when J = 1 and that 

Vj,s{u){x) = {U*Hs){x,x), xGM^ 

Putting together the previous estimates we obtain the following result. We recall that 
we have fixed a number p> A and that a = 1 — 4/p. 

Proposition 4.1 Under the previous notation and assumptions, let e = \\J — l\\i^ooi^B.^^y 
Then for any u E A and any s, < s < 1, we have, for some absolute constant C depending 
on p, 

II^J,S(^)I|C0.«(M2) < ((l + £)C/s)l|G||^l,l(R2)||u|||p(Bfl) ^ 

((l + £)C7/s)||G||^M(M2)lhlli(^(B^)- 
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The same estimate holds also for the gradient, namely 



|VPj,,H||co,.(M2) < {{l + e)C/s^)\\G\\l,2,^^2)\\u\ 



Lp{Br) 



< 



iil + e)C/s')\\G\\l,,-_ 



|2/p 



Furthermore, we have 

(4.7) W'PjA'^) - ^i,s('«)IIl°°(m2) = \\Vj-iAu)\\loo^^2) 
and, for some absolute constant C , 

(4.8) \\V{VjAu)-ViAu))\\l^ 



< \\G\\l 



< C\\G\ 



VGI 



e/s. 



Although Vj^s is nonhnear in its argument u, by a simple adaptation of the previous 
reasonings, we obtain that for any ui, U2 £ A, and for some absolute constant G depending 
on p, we have the following corresponding estimates 



\\Vj,s{ui)-Vj,s{u2)\\c^,c,^m^) < {{l + e)G/s)\\G\ 



2((l + e)C/s)||Gf^i,,(K2)i?2/P||ui 



and 



||V(T'j,,(^zi) -Pj,,(n2))||co,.(M2) < {{l + e)C/i 



\G\ 



R^/P\\u,-U2\\LriBR) < 



2iil + e)C/s')\\G\\'^2)R'^^ 



Ml - U2 



^LHBr)- 



Therefore, Vj^s '■ A — )• C^'°^{M?) is Lipschitz continuous with respect to the LP norm and 
Holder continuous with exponent 1/p with respect to the norm. 

We fix 6 such that < 6 < min{/i, ai/2}, with h, < h < 1/24, and oi > as in 
Proposition 3.7, thus depending on L only. We define the corresponding T and so as in 
We finally fix s = (A;NA)-i, < s < I/sq, and <t, < cj < s. Then we define 



Lemma 



2.2 



'Pj, 



sso y 



for any u £ A, 



where J is given by (2.3). We recall that in (2.2) we defined K = Tss^, therefore for any open 
set D C Bji we have that /^.^(xd) is the light intensity on the image plane corresponding 
to the mask D, see (2.1). 

We denote by : M — >• M the Heaviside function such that T-L^t) = for any t < and 
Hit) = 1 for any t > 0. For any constant h we set Tihit) = 1-L{t — h) for any t G M. Then, 
for any /i, < /i < 1, any s, < s < I/sq, and any a, < o" < s, we define the operator 
W : ^ ^ L°°(M2) as follows 



(4.9) 



W{u) ='Hh{IsA'^))^ for any uG^. 



Clearly, for any u £ A, W(n) is the characteristic function of an open set, which we shall 
call Vt{u). That is 

$7(ti) = {x G : Is,a{u){x) > h}, for any u e A. 



(4.10) 

In other words, Xn{u) = VV(m) = 'Hh{Is,aiu)). Moreover, whenever u 
open set contained in Bji, we shall denote i^{D) = Q{xd)- 

The final, and crucial, result of this section is the following. 



Xd, where D is an 
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Theorem 4.2 Let us fix a positive constant L. Let 5 
tion 



^2/4 and let Rq be as in Proposi- 
3.7 and such that (3.24) holds. Let us also fix R> IORq and p, p > 4:, and a = 1 — 4/p. 



We fix 6 such that < 6 < min{/i, ai/2}, with h, < h < 1/24, and ai > as in 
Proposition |3.7[ thus depending on L only. We define the corresponding T and sq as in 
We finally fix s = [kNA)^^ , < s < I/sq, and a, < a < s. Then, for any 



2.2 



Lemma 



u ^ A we define 



Is,a{u) = Vj. 



SS(-,\ 



where J is given by (2.3). 



Then for any h, 1/3 <h< 2/3, and any s, < s < I/sq, we can find positive constants 
gq, < (To < 1, and 70, depending on L, R, \\D'^{T — G)\\i^ii^^2-^, p, and ssq only, such that 
for any a, < a < aos, and any 7, < 7 < 70, the following holds. 

Let A = A^'^{r, L, R), where r = ssqRq. Let R = R + 2mi6Ro, where mi, < 
mi < 1, depends on L only. Then, for any u G A^, we have that ^(m) (s and 
r2(n) G A^'^{ri,Li,R). Here n = ssqRq < r, where Rq < 61R0/8 depends on L only, 
whereas Li > L depends on L, R, [[^^(r — G)||xi(ir2), p and ssq only. 

Moreover, the map W : A^y — )• BV{Bp^) is uniformly continuous with respect to the L^ 
norm on A^ and the distance dgt on BV{Bj^). 

Remark 4.3 We observe that the distance dgt in BV{Bj^) between W(tii) and W{u2) 
corresponds to the distance ds related to Bp^ between Q{ui) and ^(^2). 



Proof of Theorem |4.2[ The proof is a consequence of the previous analysis. We fix s, 
< s < 1/so, and h, 1/3 < h < 2/3. 

Let us begin with the following preliminary case. Let u = xd, where D ^ A, and let 
V = Tssoiu) and W = 'Hh{{Tss^{u)f). We apply Corollary 
results. 

If Ci is the open set such that W ■ 
and that {D\Bmi5r{dD)) C O and (: 



3.10 



and we obtain the following 



Xq, then, by (3.25), we notice that d^l C BmiSridD) 

^\BmMD))nn 



). Therefore Q (£ B 



R- 

We take any x £ dD and any y G M^-^i2{x), with respect to the coordinate system 
depending on x. Then we consider the points y^ = y — 5re2 and = y + 5re2. We have 
that y^ G dNs{x)\BraiSr{dD). Moreover, y~ £ D and v{y~) > 11/12, whereas y^ ^ D 
and v{y^) < 1/12. Let us call y~^ = y + to6re2, where to ^ *o depends on y, and 

v{y~^) = 1/12 whereas v{y + t6re2) < 1/12 for any t G (to, 1]. Then we use (3.26) and we 
obtain that, for any t G [— l,to]) v{y + t5re2) > 1/12 and 



^(viy + t6re2)f > 6rai/{12sso). 
at 



We may conclude that there exists a function ipi 
y = {yi,y2) e Ns{x) with \yi\ < 6ir/2, iv{y)f = 
that 



[-6ir/2,6ir/2] ] 
h if and only if y2 



such that, for any 
= ipi{yi). We recall 



(4.11) ||i;||ioo(iR2) < Ci, ||Vt;||ioo(]g2) < Ci/(sso), 
where Ci is an absolute constant, and 

(4.12) ||Vt;||co,.(M2) < C2/{ssof, 
where C2 depends on R, \\D'^{T — G)\\li(j^2\ and p only. 
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We obtain that v"^ is a C^'" function and, by the imphcit function theorem, we conclude 
that the function ipi is actually C^'". We observe that 

ll¥''lllL°°[-<5ir/2,<5ir/2] < C's, 

where C3 depends on L only. Without loss of generality, by a translation we may assume 
that (/5i(0) = 0, thus \\<fi\\L'^[-5ir/2Ar/2] < C3(5ir/2. Finally, for any ti, t2 £ [-dir/2,6ir/2], 

\ip[{ti) - ^[{t2)\ < (C4/(s5o))|tl-t2r, 

where C4 is a constant depending on L, R, ||D^(T — G)||/^i(]r2) and p only. 

Then, it is not difficult to prove that for some ri = ssqRq < r, with Rq < 61R0/8 
depending on L only, we can find Li > L, depending on L, R, \\D'^{T — G)\\ii(]^2^, p and 

sso only, such that G A^''^{ri, Li, R). Let us also remark that we have obtained that 
di{n,D) < 6r. 

Let us call e = e{s,a) = £{cr/s) = \\J — l||2,oo(B2fl)- We notice that, as a/s — )• 0"'", we 
have that e goes to as well. We also assume, without loss of generality, that e is increasing 
with respect to the variable a/s. Let us recall that, for any u E A, if w = Is,a{u), with 
< s < 1/so and < a < s, then 

(4.13) ||u;||ioo(jj2) < C5, \\Vw\\loo(^^2) < C^/isso), 
and 

(4.14) \\Vw\\co,c.(^^2) < Ce/issof, 

where C5 is an absolute constant and Cq depends on R, \\D'^{T — G)\\ii(^2-j and p only. 

For positive constants o-q, < ctq < 1, and 70, to be precised later, let us fix cr, < cr < 
o-qs and 7, < 7 < 70. We take u S A'y, v = 71<jq(u), w = Is,a{u), and D £ A such that 
11^ ~ Xd\\l^(Bb) ^ 7- Then we use Proposition 



4.1 



to infer that 



\la,s{'^) - ('7^so(Xd))^||^oo(ir2) < \\hA'^)-v'^\\L^(B.^)M\{TssQ{u)f-{Tsso{XD)f\\L^{^2) < 



||G|lii(R2)e + (2C/(sso))||G||Li(M2)||G'||vyi.i(M2)||u - XDllLP(Bfl) < (Cr/isso)) (^e + 7^^*^ 

where C is an absolute constant and consequently C7 depends on p only. 
Analogously, we can prove that 

||v (isAu) - (rs.o(xD))')|Lo.(K2) < (Cs/issof) (e + 7'/^) , 

where the constant Cg depends on \\D'^{T — G)\\ii(j^2-^ and p only. 

We now choose the positive constants fio and 70 in such a way that 

2iCr/{sso)) (e(ao) + 7o/'') <l/6 



and 



{Cs/{sso)) (e(cro) + 70^") <ai/24. 
Clearly, do and 70 depends on L, ||D^(T — G)\\Lif-^2\, p and ssq only. 
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Then we can apply to w = Is^a{u) and = i}{u) the same analysis we have used for v 
and Q in the first part of this proof. We may therefore conclude that if u S A.y and D € A 
is such that \\u — Xd||li(_b^) < 7) then Q (s Bp^, di{Q,,D) < 6r and, taken ri as before, 
possibly with a smaller Rq still depending on L only, we can find Li > L, depending on L, 
R, \\D'^{T — G)\\ii(^-^2-^, p and ssq only, such that Q G A^'°'{ri, Li, R). 

This kind of argument leads us also to show that shares the same topological properties 
of D, that is for example and have the same number of connected components of D 
and dD, respectively. 

It remains to show the uniform continuity property. We recall that the operator Vj^sso 
is Holder continuous from A, with the L^{Bji) norm, into C^'°(M2), with its usual norm. 
This means that there exists a constant C such that for any ui and U2 G A, if we call 
wi = Is,aiui) and W2 = Is,aiu2), then 

\\WI - U;2||cl.«(IR2) < C\\ui - U2\\]^^^^y 

A simple application of the previous analysis allows us to prove this claim 

Claim 1 There exists a function g : [0, +00) — )• [0,+oo), which is continuous, increasing 
and such that g{0) = 0, satisfying the following property. For any u £ A^y, for any e > 
and any x gM? we have 

(4.15) if X Beidn{u)), then \Is,a{u) - h\ > g{e). 



Let us now assume that ui and U2 belong to A--/ and let us fix e > 0. We can find 77 > 
such that if ||ni — U2\\l'^i^Bii) — ^' then ||tt;i — 'W2\\l°°{^'^) < 9{^)- 

Let us now take x G dQ{ui), that is x G such that wi{x) = h. We infer that 
1^^2(3^) — ^1 ^ ^(s)) therefore by the claim we deduce that x G -Be(9il(n2)). That is (9J7(ui) C 
B^{d^{u2))- By symmetry, we conclude that di{VL{ui), Vt{u2)) < £■ In other words, the map 
which to any u G A^./ associates the open set is uniformly continuous with r espe ct to 



3.3 



that 



the norm on A^ and the distance di. However, we have shown in Subsection 
the distances di, di, d2 and ^3 are topologically equivalent on A^'°'{ri,Li,R), to which all 
J7(u) belongs, for any u G A^y. Therefore the map A^y 9 n — )■ Vt{u) is uniformly continuous 
with respect to the norm on Ay and any of the distances di, di, d2 and ^3 related to 



We observe that 

d2(^{ui)Mu2)) = ||W(ui) - W(n2)bi(B^^^) = ||W(ni) - W{u2)\\l^b^) 

whereas 

d^{n{ui), ^{U2)) = d2{^{ui)Mn2)) + \P{^{ui)) - P{^{U2))\ = dst{W{ui),W{u2)) 

where dst is here the distance inducing strict convergence in BV{B^). Therefore we conclude 
that W : A^ — > BV{Bj0 is uniformly continuous with respect to the norm on Ay and, 
on BV{Bj0, with respect either to the norm or to the dst distance. □ 

Remark 4.4 Let us finally remark that if, instead of taking h G [1/3,2/3], we simply 
assume < h < 1, then the same analysis m ay s till be carried over. Clearly we need to 



change the values of Rq and hi in Proposition 3.7, so that they depend on h as well. As a 



consequence also the quantities introduced in the above Theorem 4.2 would depend on h. 
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5 Analysis of the inverse problem 



Throughout this section, we shall keep the notation of Theorem |4 . 2| and we shall also assume 



that the hypotheses of Theorem 4.2 are satisfied. We shall fix /i, 1/3 < h < 2/3, and s, 
< s < I/sq, and we shall take a, < a < o-qs, and 7, < 7 < 70, ao and 70 as in 
Theorem 14. 2[ 

We call the circuit to be reconstructed and we shall assume that it belongs to 



We recall that, by Proposition 3.6, A is compact with respect to the ^2 distance, which 



corresponds to the distance induced by the norm for the corresponding characteristic 



functions. Then it is an immediate consequence of the last part of Theorem 4.2, see also 



Remark 4.3 that the problem 



mmd3{n{D),no) 

admits a solution. We note that Q{D) = 0,{xd) and that here is the distance defined in 



( 3.10p related to . 

From a numerical point of view, the class A is rather difficult to handle. We try to 
reduce this difficulty by enlarging the class ^ to a class of characteristic functions of sets 
with finite perimeter. In order to keep the lower semicontinuity of the functional, we restrict 
ourselves to characteristic functions of sets with finite perimeter which are contained in A-f. 
Namely, we define the following functional Fq : A ^ [0, +00] such that for any u £ A we 
have 

(5.1) Fo{u) = dstiWiu), xno) + hP{u), 



where P is the functional defined in (3.3) with T> chosen to be Bpi^ 6 is a positive parameter 
and dst is the strict convergence distance in BV{B^. We recall that, whenever u € {0, 1} 
almost everywhere in Br and u G BV{B rj^i)^ then P{v) = P{u, Bji+i) = \Du\{Bji-^i). 
Otherwise, P{u), and consequently also Fq{u), is equal to +00. Moreover, if n G A^, in 
particular if u = xd for some D £ A, then dgtiW (u) , xno) = d3{^{u),^o), where again ^3 



is the distance defined in (|3.10|) related to Bp^ 



We look for the solution to the following minimization problem 

(5.2) m.in{Fo{u) : u G A^]. 

We notice that such a minimization problem admits a solution. 

Even if the class A^ might still be not very satisfactory to handle from a numerical point 
of view, since it somehow involves handling the class A^ we believe that from a practical 
point of view such a restriction might be dropped and we might use the class A C L^{Bji) 
instead. In fact, we have a good initial guess, given by the target circuit x^q, and it is 
reasonable to assume that the optimal mask will be a rather small perturbation of itself. 



In fact, under our assumptions, by the arguments developed in the proof of Theorem 4.2 
we can show that Vt{u) has the same topological properties of D, where xd is the element of 
A which is closest to u. Therefore if we look for a set Q.{u) as close as possible to i7o; then at 
least we need to require that the set D has the same topological properties of Oq. For this 
reason and since G A, it might be essentially the same to perform the minimization in a 
small neighbourhood of A or in the whole A. On the other hand, again by our assumptions, 
we notice that whenever the boundary of presents a corner, and this is often 
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d^lo is often the union of a finite number of segments, then rig cannot be reconstructed in 
an exact way, since Q{u), for any u G A-y, is a C^'" set, thus its boundary cannot have any 
corner. 

Besides deahng with the class Af, there are several other difficulties. In particular, 
computing -Fo(X-e) for some E C Br is not an easy task, since it involves at least the com- 
putation of the perimeters of E and of Q{xe)- Furthermore, solving a minimization problem 
in the class of sets of finite perimeter is not a straightforward task from the numerical point 
of view. 

In order to solve these difficulties, we use the following strategy. We approximate, in 
the sense of F-convergence, the functional Fq with a family of functional {-Fe}e>o which are 
easier to compute numerically and are defined on a set of smooth functions. 

As in Section 2.3, we take a C°° function : M — )• M such that (p is nondecreasing, 
(f){t) = for any t < —1/2 and (p{t) = 1 for any t > 1/2. For any r/ > and any r G M, let 



Then we have the following result. 

Proposition 5.1 For any rj > 0, let '■ A ^ C^'°'{M?) be defined as 

^n{u) = <t)n,h{Is,a{u)), for any u£ A. 

Then, for any r], < rj < h, is Holder continuous, with exponent 1/p, from A, with 
the L^{Bfi) norm, into C^'°(M^), with its usual norm. 

Furthermore, as r] ^ 0^, (W — ^ri) '■ A^ — )• BV{Bj^) converges uniformly to zero on A-y 
with respect to the distance dgt on BV{Bj^). 

Proof. The continuity property of immediately follows by the continuity of Vj^ssq and 
by the properties of (f). We just note that the Holder exponent is fixed, whereas the Holder 
constant might depend upon r/. 

About the convergence result, we begin by recalling that W(n) = 'Hh{Is,a{u)), u G A^y. 
We use Claim [l] introduced in the proof of Theorem 4.2 We call to = min(l, sup{5(s) 



s G [0,+oo)}) and sq the positive real number such that g{so) = to/2. We call : 
[0,io/2] — )• [0, So] the continuous, increasing function which is the inverse of g on such 
intervals. For any r], < rj < to, we infer that might be different from W{u){x) 



only if X G Bg-i(^^j2){dVt{u)). By estimates like (3.15) and (3.16), which are independent of 
u G ^-y, we obtain that ||(W — ^r]){u)\\ii(^B^) converges to zero, as r/ — 0"^, uniformly for 
u G A-y. 

For any t G M and any u G A^/, we call 

P{u,t) = P{{x G : /,_,(n)(x) > t},B^). 

It remains to prove that, as — )• 0+, = Jg, |V(<I>^(u))| converges to 

\D{yV{u))\{B£) = P{u,h) uniformly for u G A'y. We argue in the following way. We have 
that, for any 7], < r] < h, 

|D($,H)|(i?^)= / \V{%{u))\= [ |</>;,,(/,,.(n))||V(/,,,H)|. 
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Since (j)'^ h — ^^ ^ small enough, uniformly with respect to n G A^y, (f)'^ ^^{Is^aiu)) = 

outside Bj^, without loss of generality, we have that 



|D($,(n))|(S^)= / 0;^(I,,<,(tx))|V(/,,.(n))|. 
By the coarea formula. 



Therefore, 

VJ~oo \ V J J -1/2 

Since \D{yV{u))\{B^) = P{u,h) and J^^^j^ cf)' {s)ds = 1, we obtain that 



\Di%iumBj^)-\D{WiumBp.)\ < r'^%'{s)\P{u,srj + h) - Piu,h)\ds. 

J -1/2 



It remains to show that, as — )• 0+, sup{|P(u, t + h) — P{u,h)\ : t G [— r//2, +r//2]} 
goes to zero uniformly with respect to u G Af. Therefore the proof is concluded by using 
the following claim. 

Claim 2 There exist a positive constant 7]q and a continuous, increasing function gi : 
[0, r]o] — [0, +cxd), such that gi{0) = 0, such that for any r], < r] < rjo, and any u £ A^y, we 
have that 

snp{\Piu,t + h)-P{u,h)\ : te[-v/2,+v/2]} <gii7]). 
The proof of Claim [2] is a straightforward, although maybe lengthy to describe, conse- 



quence of the analysis developed in the proof of Theorem 4.2 We leave the details to the 
reader. We just notice that Claim [2] is a sort of generalization of Claim [T] and the arguments 
used to prove the two claims are essentially analogous. □ 

We are now in the position of describing the approximating functionals and proving the 
T-convergence result. Let us a fix a constant pi, 1 < pi < +oo, and a continuous function 
: M [0, +oo) such that W{t) = if and only if t £ {0, 1}. Let us denote by P^, e > 0, 



functional P is defined in (3.3), again with T> = Br 



the functional defined in (3.2) with p = pi, the function W and V = B^. We recall that the 



Then, for any e > 0, let us define : A ^ [0, +oo] such that for any u G A we have 
(5.3) Fe{u) = dst{%ie)iu),XnJ + bPeiu) 

where r] : [0, +oo) — t- [0, +oo) is a continuous, increasing function such that ri{0) = 0. 

By the direct method, we can prove that each of the functionals F^, e > 0, admits a 
minimum either over A or over 

The T-convergence result is the following. 
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Theorem 5.2 Let us consider the metric space {X,d) where X = A^/ and d is the metric 
induced by the L} norm. Then, as e — )• 0"^, T -converges to Fq on X with respect to the 
distance d. 

Proof. Let us fix a sequence {en}'^=i of positive numbers converging to zero as n — )• oo. 
Let, for any n G N, = -F^^. We need to prove tliat r-lim„F„ = F. Let us also remark 
that we may extend Fn and F over L^(M^) by setting them equal to +00 outside A^y. Let 
us define Pg, e > 0, and P as the functionals which are equal to the functionals and P, 



respectively, on A^y and +00 elsewhere. We recall that P^, e > 0, and P are defined in (3.2) 



and in (3.3), respectively, with p = pi and V = Br. 

We observe that, as a consequence of Proposition |5.1| and of the stability of L-convergence 
under uniformly converging continuous perturbations, it is enough to show that F-lim^ P„ = 
P, where Pn = Pe„, n G N. Let us prove this L-convergence result. 



3.1 



to construct 
such 



The L-liminf inequality is an immediate consequence of Theorem 3.1 and of the fact 
that A-y is a closed subset of V-{Br). 

For what concerns the recovery sequence, then we argue in the following way. If u G ^ is 
such that 11^* — X-dIIli(b^) < 7, for some D ^ A, then we again use Theorem 
a recovery sequence for such a function n, that is a sequence {unj^i 
that, as n — )• 00, u„ — )• n in L^{Br) and Pn{un) — )• P{u)- 

It remains to study the case when u £ dA-y and P{u) < +00. In this case, we have 
that u = xe, where E C P_r is a set of finite perimeter, and we pick D G A such that 
Wxe — Xz)||li(_Bb) = l-EADl = 7. Then at least one of these two cases must be satisfied. 
Either there exists x G P_r\-D such that 



or there exists x € D such that 



p-s>o+ \Bp{x 
\EnBp{x 



lim 



\Bp{x) 



We choose an arbitrary sequence {pj}^^^ of positive numbers such that limj pj = 0. In the 
first case, for any j G N, we choose Ej such that XEj = Xe{^ — XBp.{x))- In the second case, 
we choose Ej such that xEj = Xe{^ — XBp.{x)) + XBp.{x)- We notice that, in either cases, for 
any j G N, Ej is a set of finite perimeter such that \\xEj — Xd\\l^{Bi{) < 7- Furthermore, as 
j — )• 00 we have that xEj Xe in L^{Bpi) and P{Ej) — )• P{E), that is P{xEj) — ^ P{xe)- 
Then the proof may be concluded by following the arguments of Section 4.2 in [2] which we 



have briefly recalled in the proof of Theorem 3.1 □ 



We remark that G A, therefore we may find a family {uc}£>o such that, as e — )• 0^, 
— ^ XrJo ™ I-'^{Br) and Pe(u£) — )• P(0o)- Without loss of generality, we may assume that, 
for any e>0,0<ne<l almost everywhere in Br and that Ue G A-y. By Proposition 5.1 



we conclude that F^{Uf.) — )■ Fq{Qq) < +00. We obtain that for any eq > there exists a 
constant Ci such that 

(5.4) minPj < Ci for any e, < e < eq. 
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Obviously, the same property is shared by the minimum values of over A. 

It remains to prove that the functionals are equicoercive over A-y, that is that the 
following result holds. 

Proposition 5.3 For any eq > 0, there exists a compact subset fC of such that for any 
e, < e < Eq, we have 

min = min . 



Proof. Let us take the constant Ci as in (5.4). Let Ue G A-y, < e < eoj be such that 



Ff.{ue) = min_4 Fg. Then we observe that the set {ue}o<e<eo satisfies the properties of 



Remark 3.4 for some constant C. Therefore {ue}o<£<eo is precompact in L^[Bji) and the 



proof is concluded. □ 



Remark 5.4 With an analogous proof, the same result of Proposition 5.3 holds if we 
replace A'y with A. 



By Theorem 5.2 and Proposition 5.3, we can apply the Fundamental Theorem of F- 
convergence to conclude with the following result. 



Theorem 5.5 We have that Fq admits a minimum over A^y and 



mini*} 







lim infFe 

£-s>0+ Ay 



lim mmFgr. 

e->0+ Ay 



Let En, n gN, be a sequence of positive numbers converging to 0. For any n G N, let 
Fn = Fe„. If {un}^=i is a scqucnce contained in A'y which converges, as n ^ oo, to u £ A-y 
in L^{Bfi) and satisfies lim„F„(u„) = lim„inf_4^ Fn, then u is a minimizer for Fq on A^, 



that is u solves the minimization problem (5.2) 



We conclude with the following remark. With the notation of Theorem 5.5 , if is a 



sequence contained in A-y which satisfies lim,i Fn{un) = lim„ inf^^ Fn, then, by Remark 3.4 
we have that, up to passing to a subsequence, actually converges, as n — )• oo 

some function u E Ay in L^{Bji). 



to 



6 Discussion 



We have provided a mathematical study of the inverse problem of photolithography. The 
approach we propose is to seek an approximate solution by formulating the geometrical 
problem using a phase-field method. We further relax the hard threshold involved in image 
exposure with an approximate Heaviside function. We show that the variational problem 
for the approximate solution is well-posed. This opens a way into designing mathematically 
rigorous numerical methods. We further show that as the approximation parameter goes to 
zero, a theoretical limit, the original optimization problem involving geometry is recovered. 
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